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^ ' Abstract 

Low energy capture cross sections are calculated within a microscopic many-body approach 
using an effective Hamiltonian derived from the Argonne V18 potential. The dynamics is treated 
within Fermionic Molecular Dynamics (FMD) which uses a Gaussian wave-packet basis to represent 
the many-body states. A phase-shift equivalent effective interaction derived within the Unitary 
Correlation Operator Method (UCOM) that treats explicitly short-range central and tensor cor- 
relations is employed. As a first application the 3 He(a,7) 7 Be reaction is presented. Within the 
FMD approach the microscopic many-body wave functions of the 3/2~ and l/2~ bound states in 
7 Be as well as the many-body scattering states in the l/2 + , 3/2 + and 5/2 + channels are calculated 
as eigenstates of the same microscopic effective Hamiltonian. Finally the S-factor is calculated 
from El transition matrix elements between the many-body scattering and bound states. For 
3 He(a,7) 7 Be the S factor agrees very well, both in absolute normalization and energy dependence, 
with the recent experimental data from the Weizmann, LUNA, Seattle and ERNA experiments. 
For the 3 H(a,7) 7 Li reaction the calculated S-factor is about 15% above the data. 
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1 Introduction 

Low energy nuclear reactions play an important role in many astrophysical scenarios. If the environ- 
!h ; mental temperature is so low that the kinetic energy of the nuclei is small compared to the Coulomb 
- - ' barrier, quantum tunneling leads to reaction rates which are exponentially dropping with decreasing 
energy. Hence it becomes increasingly difficult to impossible to measure the tiny cross sections at such 
small energies and one has often to rely on extrapolations to astrophysically relevant energies. Typ- 
ically these reactions are described with potential models where the reaction partners are treated as 
point-like nuclei interacting via nucleus-nucleus potentials fitted to experimental data on bound and 
scattering states. In a microscopic ab initio picture however, the system is regarded as a many-body 
system of interacting nucleons. The many-body wave functions are fully antisymmetrized and realistic 
nucleon-nucleon interactions are used. 

In the Fermionic Molecular Dynamics (FMD) approach [H El El S] we aim at a consistent descrip- 
tion of bound states, resonances and scattering states using realistic low-momentum nucleon-nucleon 
interactions obtained in the Unitary Correlation Operator Method (UCOM) [SJ El El El [9] . Intrinsic 
many-body basis states are Slater determinants built with Gaussian wave packets as single-particle 
states. This basis contains as special cases harmonic oscillator shell model and Brink-type cluster wave 
functions. The symmetries of the system are restored by projection on parity, angular momentum and 
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total linear momentum. The many-body eigenstates of the realistic Hamiltonian are obtained in multi- 
configuration mixing calculations. FMD has already been used successfully to describe the structure 
of nuclei in the p- and sd-shell, like the Hoyle state in 12 C [ID] or halo- and cluster-structures in Neon 
isotopes [TT] . 

In this contribution we extend FMD to the continuum where many-body states have to be matched 
to phase shifted Coulomb solutions of point like nuclei. We present first results for the 3 He(a,7) 7 Be 
and 3 H(a,7) 7 Li capture cross sections using the FMD approach. The 3 He(a,7) 7 Be reaction, which is 
the onset of the 7 Be and 8 B branches of the pp II and pp III chain of hydrogen burning, has been 
studied theoretically since long either using potential models p~2j [13] or microscopic cluster models 
[IUCE] assuming 3 He+ 4 He cluster wave functions. Polarization effects were considered by including the 
6 Li+p channel in [T6l [T7] . However, consistent ab initio calculations starting from realistic interactions 
have not been possible up to now. First attempts using Variational Monte-Carlo [18] and the No-Core 
Shell Model [19] only calculated the asymptotic normalization coefficients from the bound state wave 
functions. 

2 Fermionic Molecular Dynamics 

In FMD the intrinsic many-body basis states are Slater determinants 

\Q) = A{\ qi )® ...<8> \q A )}, (1) 
using Gaussian wave packets as single-particle states 

(x\q k ) = expj- ^"^ 2 ] ® |xJ,xt> ® |&) . (2) 

The variational set Q contains the complex parameters which encode the mean positions and mo- 
menta of the wave packets, the complex width parameters a^, and the spin directions are controlled by 
Xk an d Xk- Proton and neutron are distinguished by To restore the symmetries of the Hamiltonian 
the intrinsic wave function \Q) is projected on parity, angular momentum and total linear momentum 
P = 0. 

\Q; FMK) ® \P = 0) cm = P J MK P*P P=0 \Q) (3) 




Figure 1: Density distributions of the intrinsic basis states. Top: frozen cluster configurations, bottom: 
selected polarized configurations obtained by Variation After parity and angular momentum Projection. 
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We divide the Hilbert space into an outer region where the system is represented as 3 He and 4 He 
nuclei in their ground states interacting only via the Coulomb interaction. In the inner region the 
nuclear force acts among all nucleons and the strong interactions between the approaching nuclei will 
polarize them. Therefore one has to enlarge the many-body Hilbert space beyond the microscopic 
cluster model [20] where the Hilbert space is spanned by so called "frozen" configurations that are just 
antisymmetrized products of the 4 He and 3 He clusters in their ground state placed at various distances 
R, see Fig. ffl 

In order account for polarized clusters we generate additional configurations by a variation after 
parity and angular momentum projection procedure (VAP). The energy of the projected many-body 
state is minimized with respect to the parameters of all single particle states of the intrinsic Slater 
determinant in the sub-manifold defined by a fixed mean square radius. By dialing the value of the 
radius constraint from large to small we obtain polarized clusters down to shell model like configura- 
tions. The polarized clusters are essential for the continuum scattering states and the surface region 
of the bound states. The shell model like configurations contribute mainly to bound states. Some 
polarized configurations are shown in Fig. [I] (VAP). The many-body Hilbert space is spanned by the 
non-orthogonal projected configurations given in Eq. ([3]). The task is to find a set of intrinsic states 
\Q^) that describes the physics well. Typically, we employ of the order of 50 intrinsic states which are 
taken to be the same set for bound and scattering states, of course projected on the respective J 71 ". 

The effective Hamiltonian is diagonalized in the available Hilbert space with the proper boundary 
conditions for either bound (discrete values of the energy E n ) or scattering states (E = fc 2 /(2/i)). 

# eff \rM : E n ) = E n \ J*M, E n ) (4) 

This determines the multi-configuration mixing coefficients c^- of the microscopic many-body eigen- 
states 

\J*M,E n ) = J>2 \QM;J*MK) . (5) 

aK 

The effective Hamiltonian i7 eff = T + Vqcom is obtained within the Unitary Correlation Operator 
Method (UCOM) by transforming the realistic Argonne-V18 Hamiltonian in the two-body space [5] 
This incorporates the short-range correlations into the effective interaction Kjcom and at the same time 
leaves the phase shifts of the nucleon-nucleon scattering problem and the deuteron energy unchanged. 



3 Bound and Scattering States 

The frozen states can be rewritten using RGM basis states [211 E2]. 

mr);[e± + ]£)=A [|r^®| 3 He;f)®| 4 He;0+>]^ , (6) 

where |r, £) rd denotes a state of relative motion with angular momentum £ in which the center of masses 
of the internal ground states | 3 He; | + ) and | 4 He;0 + ) are positioned at relative distance r. Internal and 
relative motion angular momenta are coupled to total J. The RGM representation can then be used 
to match the logarithmic derivative of the relative wave function of the clusters to the asymptotic 
Whittaker function for bound states or to Coulomb scattering wave functions with phase shifts: 

<*(»■);[* J*M,E) r ^ i(*i(Ar)+tanOSf (k)) G e (kr)) , k = +^E. (7) 

Using the microscopic i?-matrix approach of the Brussels group [211 E2] we solve the Schrodinger equa- 
tion with the proper boundary conditions for bound and scattering states. 
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With frozen configurations alone the 3/2~ and l/2~ states in 7 Be are bound by only 240 keV and 
10 keV, respectively. In the full Hilbert space, including the polarized configurations, we obtain binding 
energies of 1.49 MeV and 1.31 MeV with respect to the cluster threshold. In other words we reproduce 
the centroid energy but underestimate the splitting of the 3/2~ and l/2~ states. By artificially increasing 
the spin-orbit strength we observed that the total cross section depends essentially only on the centroid 
energy while the increased splitting affects only the branching slightly. The calculated 7 Be charge radius 
of 2.67 fm is in good agreement with the experimental value of 2.647(17) fm [23] . In Fig. [2] we present 
the calculated phase shifts together with the experimental data [24] . Again we see a sizable effect when 
we compare the results using only the frozen configurations with the results obtained in the full Hilbert 
space including polarization effects. 
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Figure 2: 4 He + 3 He scattering phase shifts for P-waves left and the S and D-waves right. Dashed lines 
show results using only frozen cluster configurations, solid lines show results for the full model space. 



4 Spectroscopic amplitudes 

If the many-body system has a pronounced cluster structure, as is the case in 7 Li, one can define an 
overlap function ^ (r) for the relative motion of the clusters by means of the spectroscopic amplitude 
($(r)| J* E) = ($(r); [£ | + ]^ | J^M, E) for a given J 77 and energy E as 

^(r) = / r ,2 dr'N 1/2 (r,r f ) (§{r f )\r,E) with the RGM norm kernel 7V(r,r') = ($(r)|$(r')> • (8) 




Figure 3: From left to right: overlap functions rvp^ (r) for the relative motion of the clusters for the 
3/2~ ground state and l/2 + scattering states at energies of 0.1 MeV, 1.0 MeV and 2.0 MeV. Dashed 
lines show the Whittaker and Coulomb functions matched at the channel radius a=12 fm, respectively. 
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Figure 4: Dipole strength r^/ 2 (r) • r • r^]/ 2 (r) between 3/2~ bound state and l/2 + scattering states 
with energies of 0.05, 0.1, 1.0 and 2.0 MeV. Dashed lines for Whittaker and phase shifted Coulomb 
functions matched at 12 fm. 



Fig. [3] shows three J 7r =l/2 + scattering overlap functions for different energies and the one for the 
J n =3/2~ ground state. The dashed lines are the Whittaker function and the phase shifted Coulomb 
wave functions for scattering and bound states, respectively. When there is no nuclear interaction 
between the clusters the overlap function has to go over into the Coulomb wave functions. The deviation 
from the Coulomb solution indicates that the nuclear interaction acts up to 10 fm in the scattering states 
up to about 6 fm in the more compact ground state. The ground state is an £=1 state with one node, 
while the £=0 scattering states are suppressed in the interior due to the Pauli principle and have two 
nodes. The £=0 states with zero and one node are Pauli forbidden and do not exist. 

When looking at the dipole strength displayed in Fig. 0] it can be seen that at low energies a large 
fraction of the dipole matrix element comes from the external region where only the Coulomb force acts. 
But it should also be noted that even at energies as low as 50 keV we have a significant contribution 
from the interior region, contrary to common believe that external capture is a good approximation at 
low energies. 



5 Capture Cross Section 

At low energies we can restrict ourselves to electric dipole transitions from the S- and D-wave channels. 
Using the microscopic bound and scattering wave functions and adding up all contributions from the 
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Figure 5: The astrophysical 5-factor for the 3 He(a,7) 7 Be reaction. The FMD result is given by the 
black solid line. Recent data are shown by dark (colored) symbols. Old data are included as gray 
symbols together with an i?-matrix fit. 
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J n = l/2 + , 3/2 + , 5/2 + scattering states to the J* = 3/2~, l/2~ bound states we calculate the radiative 
capture cross section which is presented in form of the astrophysical S'-factor in Fig. [51 It agrees 
very well, both in absolute normalization and energy dependence, with the recent experimental results 
obtained at the Weizmann Institute [25], the LUNA collaboration [26j, at Seattle [27] and by the ERNA 
collaboration [28] . One has to keep in mind that in this ab initio approach there are no parameters that 
can be adjusted to the data. Input are the effective realistic nucleon-nucleon interaction and a suitably 
chosen many-body Hilbert space. 

We also did the corresponding calculation for the 3 H(a,7) 7 Li capture cross sections and found that 
the measured data are about 15% lower than our calculation, see Fig. E3 
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Figure 6: The astrophysical 5-factor for the 3 H(a,7) 7 Li reaction. Recent data are shown by dark 
(colored) symbols together with an i?-matrix fit. The FMD result is given by the black solid line. Old 
data are included as gray symbols. 

Future investigations of our results have to reveal why the microscopic FMD model together with a 
realistic nuclear force describes the data while other models fail. Compared to previous works there are 
two main differences, the NN-interaction and the many-body Hilbert space. As effective interaction and 
Hilbert space are strongly correlated we do not expect to identify a single simple reason, like lacking 
momentum dependence in potential models or too simplified NN-interactions in the microscopic cluster 
models. It is however clear that on has to include polarized clusters and shell model like 7 Be ( 7 Li) states 
in order to get a successful picture. This brings about the very interesting and general question how 
composite polarizable many-body systems perform quantum tunneling. 
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